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Abstract. The first aim of the present paper is to compare various sub-Riemannian 
structures over the three dimensional sphere S 3 originating from different constructions. 
Namely, we describe the sub-Riemannian geometry of S 3 arising through its right Lie 
group action over itself, the one inherited from the natural complex structure of the 
open unit ball in C 2 and the geometry that appears when considering the Hopf map 
as a principal bundle. The main result of this comparison is that in fact those three 
structures coincide. 

In the second place, we present two bracket generating distributions for the seven 
dimensional sphere S 7 of step 2 with ranks 6 and 4. These yield to sub-Riemannian 
structures for S 7 that are not present in the literature until now. One of the distributions 
can be obtained by considering the CR geometry of S 7 inherited from the natural complex 
structure of the open unit ball in C 4 . The other one originates from the quaternionic 
analogous of the Hopf map. 

1. Introduction 

One of the main objectives of classical sub-Riemannian geometry is to study manifolds 
which are path-connected by curves which are admissible in a certain sense. In order to 
define what does admissibility mean in this context, we begin by setting notations that 
will be used throughout this paper. Let M be a smooth connected manifold of dimension 
n, together with a smooth distribution Tt C TM of rank 2 < k < n. Such vector bundles 
are often called horizontal in the literature. An absolutely continuous curve 7 : [0, 1] — >• M 
is called admissible or horizontal if j(t) £ TL a.e. 

Distributions satisfying the condition that their Lie-hull equals the whole tangent space 
of the manifold at each point play a central role in the search for horizontally path- 
connected manifolds. Such vector bundles are said to satisfy the bracket generating con- 
dition. To be more precise, define the following real vector bundles 

H 1 = H, H r+1 = [H r ,H] + U r for r > 1, 

which naturally give rise to the flag 

n = n 1 c n 2 c n 3 c . . . . 
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Then we say that a distribution is bracket generating if for all x € M there is an r(x) € Z + 
such that 



If the dimensions dim?^ do not depend on x for any r > 1, we say that TL is a regular 
distribution. The least r such that (pQ) is satisfied is called the step of TL. We will focus on 
regular distributions of step 2. In [16] the reader can find detailed definitions and broad 
discussion about terminology. 

The following classical result shows the precise relation between the notion of path- 
connectedness by means of horizontal curves and the assumption that TL is a bracket 
generating distribution. 

Theorem 1 ([91 120j). If a distribution TL C TM is bracket generating, then the set of 
points that can be connected to p € M by a horizontal path is the connected component of 
M containing p. 

Thus, the search for horizontally path-connected manifolds can be reduced to the search 
of appropriate distributions defined on them. We define the class of manifolds which will 
be of our concern. 

Definition 1. A sub-Riemannian structure over a manifold M is a pair {TL, (-, •)), where 
TL is a bracket generating distribution and (-, •) a fiber inner product defined onTL. In this 
setting, the length of an absolutely continuous horizontal curve 7 : [0, 1] — ► M is 



Jo 

where ||7(t)|| 2 = ( A f(t),j(t)) whenever j(t) exists. The triple (M,TL, (•,•)) is called sub- 
Riemannian manifold. 

Thereby, restricting ourselves to connected sub-Riemannian manifolds endowed with 
bracket generating distributions, it is possible to define the notion of sub-Riemannian 
distance between two points. 

Definition 2. The sub-Riemannian distance d(p, q) € M>o between two points p,q G M 
is given by d(p,q) := mfl^), where the infimum is taken over all absolutely continuous 
horizontal curves joining p to q. 

An absolutely continuous horizontal curve that realizes the distance between two points 
is called a horizontal length minimizer. 

Remark: The connectedness of M and the fact that TL is bracket generating, assure that 
d(p, q) is a finite nonnegative number. Nevertheless, the bracket generating hypothesis, 
required for the previous definition, is possible to be weakened. In fact, in [22J the au- 
thor finds a necessary and sufficient requirement to horizontal path-connectedness for a 
manifold in terms of the corresponding distribution. Clearly, this theorem contains, as a 
particular case, the bracket generating condition. 



(1) 



TL r x {x) = T X M. 




Historically, the first examples of sub-Riemannian manifolds that have been considered 
were Lie groups, see e.g. [2J [5j 13 HSJ E], because due to its algebraic structure, it is 
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possible to reduce the problem of finding globally denned bracket generating distributions 
to finding such distributions at the identity by considering right (or left) invariant vector 
fields. An important role has been played by considering domains in Euclidean spaces 
with special algebraic structures (such as the Heisenberg groups, H— type groups as their 
natural generalizations to Clifford algebras, Engel groups, Carnot groups, etc.). Particular 
attention have had the three dimensional unimodular Lie groups which were studied, for 
example, in [2], [5] and [12]. The main purpose of this communication is to present recent 
results concerning different sub-Riemannian structures of the second simplest family of 
examples of manifolds, namely, spheres. We also be inspired by the article [23] . where the 
close relation between the Hopf map and physical applications is presented. 

The following celebrated theorem in topology, see [I], gives a very strong restriction on 
the problem of finding globally defined sub-Riemannian structures over spheres. 

Theorem 2 (|lj). Let S n ~ l = {x G M n : ||x|| 2 = 1} be the unit sphere in M. n , with respect 
to the usual Euclidean norm \\ ■ ||. Then S" 1-1 has precisely g(n) — 1 globally defined and 
non vanishing vector fields, where g(n) is defined in the following way: if n = (2a + l)2 b 
and b = c + Ad where < c < 3, then g(n) = 2 C + 8d. 

In particular, two classical consequences follow: S , S 3 and S 7 are the only spheres 
with maximal number of linearly independent globally defined non vanishing vector fields, 
and the even dimensional spheres have no globally defined and non vanishing vector fields. 

The condition that a manifold M has maximal number of linearly independent globally 
defined non vanishing vector fields is usually rephrased as saying that M is parallelizable. 
The fact that S l , S 3 and S 7 are the only parallelizable spheres was proven in [6] and that 
the even dimensional spheres have no globally defined and non vanishing vector fields is a 
consequence of the Hopf index theorem, see |24j . 

This theorem permits to conclude at least two major points of discussion: there is no 
possible global basis of a sub-Riemannian structure for spheres with even dimension and 
it is impossible to find a globally defined basis for bracket generating distributions, except 
for S 3 and S 7 . The fact that S 3 and S 7 can be seen as the set of quaternions and octonions 
of unit length will play a core role in many arguments throughout this paper. 

The main results that we present here are: a comparison between three sub-Riemannian 
structures of S 3 , namely, the one arising through its right Lie group action over itself as 
the set of unit quaternions, the one inherited from the natural complex structure of the 
open unit ball in C 2 and the geometry that appears when considering the Hopf map as 
a principal bundle. Notice that this structure admits a tangent cone isomorphic to the 
one dimensional Heisenberg group in a sense of Gromov-Margulis-Mitchell-Mostow con- 
struction of the tangent cone [131 118} [19]. Considering CR-structure of S , inherited from 
the natural complex structure of the open unit ball in C 4 , we obtain a 2-step bracket 
generating distribution of the rank 6. This construction intimately related to the Hopf 
fibration S 1 — ► S 7 — > CP 3 . Making use of the quaternionic analogue of the Hopf map 
S 3 — > S 7 — > S A , we present another 2-step bracket generating distribution that has the 
rank 4. We conclude that the sphere S 7 admits two principally different sub-Riemannian 
structures. The tangent cone, in the first case, is isomorphic to the 3-dimensional Heisen- 
berg group, and in the second case it has a structure of the quaternionic Heisenberg-type 
group with 3-dimensional center [;7j. In both cases we present the basis of the horizontal 
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distribution that is very useful in future studies of geodesies and hypoelliptic operators 
related to the spherical sub-Riemannian manifolds. 

2. S 3 AS A SUB-RlEMANNIAN MANIFOLD 

Throughout this paper H will denote the quaternions, that is, M = (M 4 , +, o) where + 
stands for the usual coordinate-wise addition in M 4 and o is a non-commutative product 
given by the formula 

(x + x\i + x 2 j + x 3 k) o (y + yi i + y 2 j + y 3 k) = 
= (xoVo ~ x\yi - x 2 y 2 - x 3 y 3 ) + (xiy + x yi - x 3 y 2 + x 2 y 3 )i+ 
+(x 2 yo + x 3 yi + x$y 2 - xiy 3 )j + (x 3 y - x 2 yi + x x y 2 + x y 3 )k. 

It is important to recall that H is a non-commutative, associative and normed real 
division algebra. Let q = t + ai + bj + ck £ M, then the conjugate of q, is given by 

q = t — ai — bj — ck. 

We define the norm \q\ of q 6 H, as for the complex numbers, by \q\ 2 = qq. 

Given the fact that the sphere S 3 can be realized as the set of unit quaternions, it has 
a non abelian Lie group structure induced by quaternion multiplication. 

The multiplication rule in EI induces a right translation R y {x) of an element x = xq + 
x\i + x 2 j + x 3 k by the element y = yo + y\i + y 2 j + y 3 k. The right invariant basis vector 
fields are defined as Y(y) = Y(0)(R y (x))*, where Y(0) are the basis vectors at the unity 
of the group. The matrix corresponding to the tangent map (Ry(x))*, obtained by the 
multiplication rule, becomes 



{R y (x))i 



( yo yi yi 2/3 \ 

-yi yo -y3 yi 

-2/2 yz yo -yi 

V -2/3 -2/2 2/1 2/0 / 



Calculating the action of (R y {x)) if in the basis of unit vectors of M 4 we get the four 
vector fields 

N{y) = yod yo + yid yi + y 2 d V2 + y 3 d V3 
V{y) = -yid yo + y d yi - y 3 d y2 + y 2 d y3 
X{y) = -y 2 d yo + y 3 d yi + y d m - yxd yz 
Y{y) = -y 3 d yo - y 2 d yi + y x d m + y d y3 . 

It is easy to see that N(y) is the unit normal to S 3 at y € S 3 with respect to the usual 
Riemannian structure (•, •) in TK 4 . Moreover, for any y G S 3 

(N(y), V(y)) y = (N(y), X(y)) y = (N(y), Y(y)) y = 

and 

(N(y),N(y)}y = (V(y),V(y)) y = (X(y),X(y)) y = (Y(y),Y(y)) y = 1. 
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Since the matrix 

(-2/1 2/0 -2/3 2/2 \ 
-2/2 2/3 2/0 -2/1 
-2/3 -2/2 2/1 2/0 / 
has rank three, we conclude that the vector fields {V(y), X(y),Y (y)} form an orthonormal 
basis of T y S 3 with respect to (-, -) y . 

Observing that [X, Y] = 2V, we see that the distribution span{X, Y} is bracket gener- 
ating, therefore it satisfies the hypothesis of Theorem [TJ The geodesies of the left invariant 
sub-Riemannian structure of S 3 are determined in [8], while in [15] the same results are 
achieved by considering the right invariant structure of S 3 . 

Notice that the distribution span{X, Y} can also be defined as the kernel of the contact 
one form 

uj = -yi dy + y dy% - y 3 dy 2 + 2/2 dy 3 . 

Remark: It is easy to see that [V, Y] = 2X and [X, V] = 2Y, therefore the distributions 
span{y, V} and span{X, V} are also bracket generating. The corresponding contact forms 
are 

6 = -y 2 dy + y 3 dyi + y dy 2 - y\ dy 3 

and 

V = -2/3 dy - y 2 dy 1 + y 1 dy 2 + y dy 3 

respectively. This means that there is a priori no natural choice of a sub-Riemannian 
structure on S 3 generated by the Lie group action of multiplication of quaternions. Any 
choice that can be made, will produce essentially the same geometry. 

3. S 3 AS A CR MANIFOLD 

Consider S 3 as the boundary of the unit ball B A on C 2 , that is, as the hypersurface 

S 3 := {(z, w) € C 2 : zz + ww = 1}. 

The sphere S 3 cannot be endowed with a complex structure, but nevertheless it possess 
a differentiable structure compatible with the natural complex structure of the ball B as 
an open set in C 2 . We will show that this differentiable structure over the sphere S 3 (CR 
structure) is equivalent to the sub-Riemannian one considered in the previous section. We 
begin by recalling the definition of a CR structure, according to [4j. 

Definition 3. Let W be a real vector space. A linear map J : W — > W is called an almost 
complex structure map if J o J = —I, where I : W — > W is the identity map. 

In the case W = T p M 2n , p = (xi,yi, . . . , x n , y n ) G K 2n , we say that the standard almost 
complex structure for W is defined by setting 

J n (d Xj ) = d yv Jn{dy 3 ) = -d Xj , l<j<n. 

For a smooth real submanifold M of C n and a point p £ M, in general the tangent 
space T p M is not invariant under the almost complex structure map J n for T p (C n ). We 
are interested in the largest subspace invariant under the action of J n . 
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Definition 4. For a point p £ M , the complex or holomorphic tangent space of M at p 
is the vector space 

H p M = T P M n J n (T p M). 
In this setting, the following result takes place (see [1]). 
Lemma 1. Let M be a real submanifold ofC n of real dimension 2n — d. Then 

2n-2d< diniR H p M < 2n - d, 
and diniK H p M is an even number. 

A real submanifold M of C n is said to have a CR structure if diniR H p M does not 
depend on p £ M. In particular, by the previous lemma, every smooth real hypersurface 
S embedded in C n satisfies dim^ M p S = 2n — 2, therefore S is a CR manifold. This fact 
is applied to every odd dimensional sphere. 

The question addressed now is to describe H p S 3 . By the discussion in the previous 
paragraph, H P S 3 is a complex vector space of complex dimension one. This description 
can be achieved by considering the differential form 

to = zdz + wdw 

and observing that kerw is precisely the set we are looking for. Straightforward calcula- 
tions show that 

kerw = sp&n{wd z — zd w }. 
In real coordinates this corresponds to 

wd z - zd w = -A-X + iY), 

where X and Y were defined in Section [21 It is important to remark that this is precisely 
the maximal invariant J2— subspace of T p S 3 , namely 

J 2 (X)=Y, J 2 (Y) = -X, 

then J 2 (span{X,y}) = span{X,y}, but J 2 {V) = -N $ T P S 3 for any point p G S 3 . 
Therefore, the distribution corresponding to the right invariant action of S 3 over itself is 
the same to its one dimensional (complex) CR structure. 

Remark: The distribution associated to the anti-holomorphic form 

lj = zdz + wdw 
is the same as the previous one. More explicitly: 

kerw = span{— wdz + zd^} 
and in real coordinates this corresponds to 

-wdg + zde = ^(X + iY). 

The same almost complex structure as the previously described can be obtained by 
means of the covariant derivative of S 3 considered as a smooth Riemannian manifold 
embedded in IR 4 . Namely, in [15] it is introduced the mapping J(X) = DxV, were D 
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denotes the Levi-Civita connection on tangent bundle to S 3 and V is the vector field 
defined in Section [2l 

4. S 3 AS PRINCIPAL BUNDLE 

In this section we describe how the structure of a principal S 1 — bundle over S 3 induces a 
bracket generating distribution on S 3 . Namely, it is possible to consider S 3 as a S 1 — space, 
according to the action 

A • (z, w) = (Xz, Xw), 

where A G S 1 = {v G C : \v\ 2 = 1} and (z, w) G S 3 = {{z, w) G C 2 : \z\ 2 + \w\ 2 = 1}. 
Consider the Hopf map h : S 3 — > S 2 as a principal S 4 — bundle, see [Hj, given explicitly 

by 

h(z,w) = (\z\ 2 — \w\ 2 , 2zw), 

where S 2 = {(x, Q G R x C : x 2 + |C| 2 = 1}- Clearly, fa is a submersion of S 3 onto S 12 , and 
it is a bijection between S 3 /S 1 and S 2 , where S 3 /S 1 is understood as the orbit space of 
the S 1 — action over S 3 , previously defined. 

Let p = (xq, Co) ^ S 2 . It is easy to verify that h _1 {p) = (zq, ^o) mod S , where (zq, wq) 
is one preimage of p under h. Consider the great circle 

lp (t) = e 2 ™ t (z ,w ), te[0,l], 

in S 3 , that projects to p under the Hopf map. Now consider the following vector field 
associated to j p , defined by the tangent vectors 

%(t) = 2TTie 2 ™\z ,w ) e T lp(t) S 3 . 

We write the Hopf map in real coordinates, where, in particular, 7 p (i) = (z(t),w(t)) = 
(xo(t) + ixi(t), X2(t) + ixs(t)) = (xo(t),xi(t),X2(t),xs(t)). It is easy to see that 

/ x (t) xi(t) -x 2 {t) -x 3 (t) 

i d J P (t) h } = 2 I x 2(t) X 3 (t) X Q (t) Xl(t) 

\ -x 3 (t) x 2 (t) Xl(t) -x (t) 
Thus, the Hopf map induces the following action over the vector field j p (t): 

I x (t) Xl(t) -x 2 (t) 
\ d -y P {t)h]i P {t) = 4vr x 2 (t) x 3 (t) x (t) 

\-X 3 (t) X 2 (t) Xl(t) 

Therefore, if [d 7 ^h] is a full rank matrix, we would have characterized the kernel of it, 
by 

(2) kerd 7p(t) /i = span{7 p (t)}. 

Notice that, using the notation of Section [21 the following identity holds 

(3) %(t) = 2irV(7„(t)). 

There is a simple way to see that the matrix [d 7p (t)/i] is full rank. Denote by Di the 3x3 
matrix obtained by deleting the i— th column of [dj p r t \h], i = 1, 2, 3 or 4, then 

det(£>i) 2 + det^a) 2 + det(D 3 ) 2 + det(D 4 ) 2 = 








( \ 


x x (t) 







x 2 (t) 
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= ((x (t)f + { Xl (t)f + (z 2 (i)) 2 + (x 3 (t)ff = 1 

implies that [d^Mh] is full rank. 

Before describing how the Hopf map induces a horizontal distribution, it is necessary 
to present some definitions. 

Definition 5 (Ehresmann Connection [19] . chapter 11). Let M and Q be two differ -entiable 
manifolds, and let ir : Q — > M be a submersion. Denoting by Q m = 7r _1 (m) the fiber 
through m G M , the vertical space at q is the tangent space at the fiber Q^^q) and it is 
denoted by V q . 

An Ehresmann connection for the submersion ir : Q — ► M is a distribution TL C TQ 
which is everywhere transversal to the vertical, that is: 

V q ®TL q = T q Q. 

We apply Definition [5] to S 3 in order to define the Ehresmann connection. Since we 
know that kevd p h = sp&n{V(p)}, for every p € S 3 by ([2]) and ([3j), and moreover, 

{X(p),V(p)) p = (Y(p),V(p)) p = {X(p),Y(p)) p = 

where (•, -) p stands for the usual Riemannian structure defined at p G S 3 , we see that 

(4) H p = s V( m{X(p),Y(p)} 

is an Ehresmann connection for the submersion h : S 3 — > S 2 with V(p) as a vertical space. 

Definition 6. Let G be a Lie group acting on Q and ir : Q — > M a submersion, with 
Ehresmann connection ri, which is a fiber bundle with fiber G. The submersion tt is called 
a principal G— bundle with connection, if the following conditions hold: 

• G acts freely and transitively, 

• the group orbits are the fibers of ir : Q — > M (thus M is isomorphic to Q/G and it 
is the canonical projection) and 

• the G— action on Q preserves the horizontal distribution Tt; 

We conclude that the Hopf fibration is a principal S 1 — bundle with connection TC, defined 
in gl). 

Definition 7. A sub-Riemannian metric {TL, (•, •)) on the principal G-bundle tt : G — > M 
is called a metric of bundle type if the inner product (•, •) on the horizontal distribution TL 
is induced from a Riemannian metric on M. 

The sub-Riemannian metric (•, -)\hi obtained by restricting the usual Riemannian metric 
of S 3 to the distribution TL is, by construction, a metric of bundle type. 

Thus, the Hopf map indicates in a very natural topological way how to make a natural 
choice of the horizontal distribution TL that was not obvious when we considered the right 
action of S 3 over itself. 

Remark: Observe that the considered vector fields coincide with the right invariant vector 
fields. This phenomenon does not appear when we change the right action to the left action 
of S 3 over itself. 
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5. Tangent vector fields for S 7 

In Sections 5 - 7 we study different sub-Riemannian structures over the sphere S 7 , 
making use of ideas of Sections 2-4. As a result, we obtain two principally different types 
of horizontal distributions. One of them of rank 6 and other of the rank 4. Moreover, as 
we shall see, the sub-Riemannian structure induced by the CR-structure and quaternionic 
analogue of the Hopf map are essentially different. We start from the construction of the 
tangent vector fields to S 7 . 

The multiplication of unit octonions is not associative, therefore S 7 is not a group in 
a contrast with S 3 . Nevertheless, we still able to use the multiplication law in order to 
fugue out global tangent vector fields. To do this, we present a multiplication table for the 
basis vectors of M 8 . The non-associative multiplication gives rise to the division algebra 
of the octonions 

O = span{e ,ei,e2,e3,e4,e5,e 6 ,e 7 }. 





eo 


ei 


e2 


e3 


e 4 


es 


ee 


e7 


eo 


eo 


ei 


e2 


e3 


e 4 


es 


ee 


e7 


ei 


ei 


-eo 


e3 


-e 2 


es 


-e 4 


-e 7 


ee 


e2 


e2 


-e3 


-eo 


ei 


ee 


e7 


-e 4 


-es 


e3 


e3 


e2 


-ei 


-eo 


e7 


-ee 


es 


-e 4 


e 4 


e 4 


-es 


-ee 


-e 7 


-eo 


ei 


e2 


e3 


es 


es 


e 4 


-e 7 


ee 


-ei 


-eo 


-e3 


e2 


ee 


ee 


e7 


e 4 


-es 


-e2 


e3 


-eo 


-ei 


e7 


e7 


-ee 


es 


e 4 


-e3 


-e 2 


ei 


-eo 



Table 1. Multiplication table for the basis of O. 



According to the Table HJ the formula for the product of two octonions is presented in 
the Appendix, subsection 18.11 This multiplication rule induces a matrix representation of 
the right action of octonion multiplication, given explicitly by: 



2/0 


-2/1 


-2/2 


-2/3 


-2/4 


-2/5 


-2/6 


-2/7 


2/1 


2/0 


2/3 


-2/2 


2/5 


-2/4 


-2/7 


2/6 


2/2 


-2/3 


2/0 


2/1 


ye 


2/7 


-2/4 


-2/5 


2/3 


2/2 


-2/1 


2/0 


2/7 


-2/6 


2/5 


-2/4 


2/4 


-2/5 


-2/6 


-2/7 


2/0 


2/1 


2/2 


2/3 


2/5 


2/4 


-2/7 


2/6 


-2/1 


2/0 


-2/3 


2/2 


2/6 


2/7 


2/4 


-2/5 


-2/2 


2/3 


2/0 


-2/1 


2/7 


-2/6 


2/5 


2/4 


-2/3 


-2/2 


2/1 


i/o 



We are able to find globally defined tangent vector fields which are invariant under 
this action. We proceed by the analogy with Section 2. The explicit formulae are given 
in subsection 18.21 of the Appendix. The vector fields {Yq, . . . ,Yj} form a frame. More 
explicitly, we have that the following identity holds 

(Y i (y),Y j (y)) y = S ij , y G S 7 , i, j G {0, 1, . . . , 7}, 
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where (■, ■} is the standard Riemannian structure over R 8 . 

6. C.R- STRUCTURE AND THE HOPF MAP ON S 7 

In the book [19] it is briefly discussed the general idea of studying a sub-Riemannian 
geometry for odd dimensional spheres via the higher Hopf fibrations. Namely, consider 

S 2n+1 = | z G C n+1 . || z ||2 = ^ then ^ S<l_ act i orl on ^n+l giyen by 

A - (z ,...,z n ) = (Xz ,...,Xz n ), 
for A € S 1 and (zq, . . . ,z n ) € S 2n+ , induces the well-known principal S 1 — bundle S 1 — » 

5 2n+l C pn given explicitly by 

5 2n+1 9(z ,...,z n )M[z :...: 4 ]eCP n 

where [20 : • • • : z n] denote homogeneous coordinates. This map is called higher Hopf 
fibration. The kernel of the map h: S 2n+1 — > CP n produces the vertical space and a 
transversal to the vertical space distribution gives the Ehresmann connection. We show 
that the vertical space is always given by an action of standard almost complex structure 
on the normal vector field to S 2n+ , and the Ehresmann connection coincides with the 
holomorphic tangent space at each point of S 2n+1 . 

Theorem [2] asserts that any odd dimensional sphere has at least one globally defined 
non vanishing tangent vector field. If the dimension of the sphere is of the form 4n + 1, 
then it has only one globally defined non vanishing tangent vector field. In the case that 
the dimension of the sphere is of the form 4n + 3, then the sphere admits more than one 
vector field. The sphere S 2n+1 possess the vector field 

V n +i(y) = -yidy + y d yi - y 3 d y2 + y 2n +2dy 2n+1 + y2n+\d y2n+2 . 

Observe that this vector field has appeared already in two opportunities: the vector field 
V in Sections El [3] and H] corresponds to V2; and the vector field Y\ in Subsection 18.21 of 
the Appendix corresponds to V4. 

The vector field V n +i has the remarkable property that it encloses valuable informa- 
tion concerning the CR structure of S 2n+1 . We know by Lemma Q] that, as a smooth 
hypersurface in C n+1 the sphere S 2n+l admits a holomorphic tangent space of dimension 

dim R H p S 2n+1 = 2n 

for any point p G 5* 2n+1 . The following lemma implies the description of H p S 2n+l as the 
orthogonal complement to V^+i. 

Lemma 2. Let W be an Euclidean space of dimension n + 2 ; n > 1, and inner product 
(■,-)w- Consider an orthogonal decomposition W = span{X, Y} ©j_ W with respect to 
(•, -)w and an orthogonal endomorphism A : W — > W such that 

^(span{X,y}) = span{X,y}, 

then W is an invariant space under the action of A, i.e. 

A(W) = W. 
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Proof. Let v £ W, then for any a, (3 S K it is clear that 

(Av,aX + (3Y) W = (v,A\aX + (3Y)) W = (v, A' 1 (aX + f3Y)) w . 

Since ^4(span{X, Y}) = span{X, Y}, there exist real numbers a, b such that A~ l (aX + 
j3Y) = aX + 6Y, and therefore 

(Av, aX + (3Y) W = (v, aX + bY) w = 

which implies that Av € VK. □ 

As an application of Lemma [21 it is possible to obtain an explicit characterization of 
H p S 2n+1 . 

Lemma 3. For any p € S 2n+l , the vector space H p S 2n+1 is the orthogonal complement 
to the vector V n+ \{jp) in T p S 2n+l . 

Proof. Consider the vector space 

W p = span{iV n+1 (»} e± T p S 2n+1 * T p R 2n+2 , 

where N n+ \ (p) is the normal vector to S 2n+1 at p. The standard almost complex structure 
map J n +i : W p — ► W p is clearly orthogonal. Moreover 

J n +i{V n+ i{p)) = -N n+1 (p), J n+1 (N n+1 (p)) = V n+ i(p). 

Using the decomposition W p = W p ©j_ span{T4, + i(p), N n+ i(p)} it is possible to apply 
Lemma[2]in order to conclude that W p , which is the orthogonal complement to V n +i(p) in 
T p S 2n+1 , is invariant under J„+i. Since dimjR W p = 2n, we conclude that W p = H p S 2n+1 . 

□ 

Remark: The space HS 2n+1 can also be described as the kernel of one form 

uj = zqcIzq + . . . + z n dz n . 
Indeed, take X € HS 2n+1 , then by straightforward calculations we have 

u{X) = (X, N n+1 ) + i(X, K+i) = 0. 

Lemma [3] provides a horizontal distribution of rank 2n for the spheres S 2n+1 , by con- 
sidering the holomorphic tangent space. The question now is whether this distribution 
is bracket generating. The bracket generating condition for S 3 is already discussed in 
Section [2j Here we present a basis for bracket generating distribution of rank 6 for S 7 . 
The following ideas were introduced to the authors by Prof. K. Furutani in a private 
communication. 

Theorem 3. The subbundle TL = spanjl^, • • • >^7} ofTS 7 is a bracket generating distri- 
bution of rank 6 and step 2. 

Proof. Define the following vector fields 

v 4 i(y) = -yidy + y 5 d yi + y d V4 - y±d y5 , 

V42(y) = yedy 2 - y 7 d y3 - y 2 d y(j + y 3 d y7 , 

V5i(y) = -ynd yo - yid yi + yid yA + y d y5 , 

V52{y) = -yid m - y e d m + y^d m + y d yr , 
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and observe that V41 + v^ 2 = Y4 and U51 + v^ 2 = Y5. By straightforward calculations we 
see that 

(v 4 i (y),Y (y)) y = (v 4 2(y),Y (y)) y = (v 51 (y),Y (y)) y = (v 52 (y),Y (y)) y = 

(«4i(y),H(y))„ = («42(y), n(i/)>„ = (i; 5 i(y),n(y)) y = (^ 2 (y),^i(y)) 2/ = 

which implies that U41, U42, U51, V52 £ span{Y2, • • • ,Y-j}. The following commutation rela- 
tion 

[l>41,l>5l] + [^42,^52] = -2Yi 
implies that the distribution TL is bracket generating of step 2. □ 

Remark: It is possible to repeat the previous argument with other pairs of vector fields. 
For example, if instead of I4 and Y§ we employ Y 2 and I3, we can consider the vector 
fields 

V2i(y) = -yidy + y 5 d yi + y d y4 - yid y5 , 
V22(y) = yad y2 - y 7 d ys - y 2 d ye> + y 3 d V7 , 
V3i(y) = -y^dyo - y\d yx + y x d m + y d y5 , 
V32(y) = -yid y2 - yedy 3 + ysd y6 + y d V7 . 
We can proceed in a similar way if we use Yq and Y7. 

We conclude this section by proving that the line bundle span{K„ + i} is the vertical 

space for the submersion given by the Hopf fibration S 1 <—* S 7 — CP". This implies 
that the distribution TL defined in Theorem [3] is an Ehresmann connection for h. To 
achieve this, we recall that the charts defining the holomorphic structure of CP n are given 
by the open sets 

U k = {[z : ■ ■ ■ : z n ] : z k ^ 0}, 
together with the homeomorphisms 

<Pk ■ U k -» C" 

Then, without loss of generality we will assume that n = 3 and we will develop the explicit 
calculations for k = 0. The other cases can be treated similarly. 
Using the chart (Uq, cpo) defined above, we have the map 



tpo o h : S 
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(z , Zl ,z 2 ,z 3 ) » 



which in real coordinates can be written as 

,„„ „ h(<r n rr-~\ — I XQX2+XXX3 X0X3-X1X2 X0X4+X1X5 X0X5-X1X4 XQXQ+XIX7 XQX7-X1XS 

(f0 O U[Xq, . . . , Xj) — I — 2,2 — , — T 2 , r 2 j — T 2 , T 2 ) — „2 , _2 > — T 2 , '2 — j — _2 , T 2 — 

\ 1 01 1 ' 1 ' 1 01 

The differential of this mapping is given by the matrix 



d((fo oh) = 
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/ 


(xl 




(*§ 




XQ 


xi 














(x? 






(xg+xf)^ 


X Q ~\~ X 


X Q ~\~ X 




-Xq)x 3 -\-2x xix 2 




-Xq)x 2 —2x xix 3 


XI 


XQ 






















(x 2 +x-(P 


x Q -{-x 1 


x Q +x 1 




(x? 


-Xq)x4 — 2a;Qa;iX5 




— x\ )x§ —2xo x 1 x 4 








x a 


xi 












(xl+x{) 2 




(*g + *?) 2 


X Q ~\~ X 


X Q ~\~ X 




(x? 


-X^X^^XQXlX^ 




-Xq)x4 — 2xqxiX5 








xi 


XQ 












(xl + xl) 2 




(x^+xir 


X q-\-X 


X q~\- X ^ 




(x? 


~x'^ ) )xq~2xqx\X7 


(*n 


-x\)xj — 2xqx\xq 














x 


x 1 






{xl+xD* 




(x%+xjr 


^0 1 


Xq ~\~X 


V 




-x 2 ,)x 7 + 2x x 1 x 6 




-Xq)xq — 2xqx\x^ 














xi 


XQ 




{xl+xD* 




(x%+xir 


Xq ~\~x ^ 


X Q~\~ X 



By straightforward calculations, we know that 

det([d(p o H)} [d(<p o H)Y) = (x 2 + xjy 8 , 
therefore, the matrix d(tpo ° H) has rank 6 and the dimension of its kernel is 2: 

diniK ker d(ipo o H) = 2. 

Moreover, since 

d(p o #)(iV„ +1 ) = d(<A) o F)(K+i) = 0, 
by direct calculations, we conclude 

ker d(if o H) = span{AT„ +1 , V n+ i}. 

This implies that 

ker dH = span{V^+i}. 

7. Application of the first quaternionic Hopf fibration 

Trying to imitate the work already done for S 3 , we find through the quaternionic Hopf 
bundle S" 3 — > S —> S 4 one of the natural choice of horizontal distributions. We consider 
the quaternionic Hopf map given by 



(5) 



,2zw) ' 



x§ 



r 2 

Xq 



h : S 7 -► S 4 
(z, w) i — ► (\z\ 2 — \w[ 

It can be written in real coordinates: 

(6) h(xo, ■ ■ . , x 7 ) = (xq + x\ + x\ + x\ - x\ - xg - x^ - x\ 

2(x X 4 + X1X5 + X 2 X 6 + X 3 X 7 ), 2(-X X 5 + X1X4 - X 2 X 7 + X 3 X 6 ), 

2(-xox 6 + xix 7 + X2X4 - X3X5), 2(-xqX7 - xixe + X2X5 + X3X4)). 
The differential map dh of h is the following: 

/ 



dh 



\ 



x 


Xl 


^2 


^3 


— X4 


-X5 


X4 


X5 


X 6 


X 7 


x 


Xl 




X4 


-x 7 


x 6 


Xl 


-Xq 


-Xq 


X 7 


X4 


-X5 


X2 


-X3 


-x 7 


-x 6 


x 5 


X4 


X3 


X2 



-Xq 
X 2 
X3 
-Xq 
-Xl 



-X 7 \ 

X3 
-X2 
Xl 

-XQ ) 



Since no one of the commutators [Y$,Y^], i,j = 1,...,7 coincides with the k = 
1, . . . , 7, we look for the kernel of dh among the commutators Yy, i,j = 1, . . . , 7. We 
found that [d/i]Y 4 5 = [^1*46 = [d/i]Y 5 6 = 0. Define V = {I45, Y 46 , Y 56 }. 
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Our next step is to find the horizontal distribution span{7Y} that is transverse to 
span-jV} and bracket generating: span{7^} p © span{V} p = T P S 7 for all p € S 7 . To begin 
we define five basis for horizontal distributions, that we shall work with. The numeration 
is valid only for this section. 

H = {Y 47 ,Y 57 ,Y 67 ,W}, 

Hi = {y 34 , y 35 , y 36 , 137}, n 2 = {Y 24 , Y 25 , y 26 , Y 27 }, 

H 3 = {Y u , Y 15 , Fie, Y 17 }, H 4 = {Y 04 , Y 05 , Y m ,Y 07 }, 

where the vector field W will be defined later and the notation Yofc = Yj, is chosen for the 
convenience. 

We collect some useful information about sets H m , m = 0, . . . , 4, that we exploit later. 

1. All vector fields inside H m , m = 0,1,2,3,4 are orthonormal (we do not count W 
before we precise it). 

2. All of collections Ti m , m = 0,1,2,3,4 are bracket generating with the following 
commutator relations: 

^[Y j4 , Y j5 ] = y 45 , \[Yj4, Y j6 ] = y 46 , ^[Y j5 , Y j6 ] = y 56 , j = o, 1, 2, 3, 

\[Yai, y 57 ] = Y 45 , ^[Y 47 , Y 67 ] = Y m , l\y 57 , Y %7 ] = y 56 - 

3. We aim to calculate the angles between the vector fields from T~L m , m = 0, 1, 2, 3, 4 
and between vector fields from TC m and V. Beforehand, we introduce the following 
notations for the coordinates on the sphere S 4 given by the Hopf map S 3 — > S 7 — > 

aoo = vl + yl + vl + yl-yl-yl-vl- y 7 , 

an = 2(y y4 + 2/1^5 + ^2^6 + 2/32/7), 
(7) a 22 = 2(-yQy<s + yiy 4 - y 2 y 7 + y 3 y & ), 

033 = 2(-y 2/6 + 2/12/7 + 2/22/4 - 2/32/5), 
a 44 = 2(-y 2/7 - 2/12/6 + 2/22/5 + 2/32/4) • 

The first index of a m k reflects the number of the collection 7i m , where they will 
appear and the second one is related to the number of the coordinate on S 4 . 

We start from Tto and calculate the inner products: 

0^45, 5-67) = ~(Y 46 ,Y 57 ) = (Y 56 ,Y 47 ) = a 00 . 
All other vector fields are orthogonal. We continue for Hi. 

(Y 45 ,Y 36 ) = -(Y 46 ,Y 35 ) = (Y 56 ,Y 34 ) = an 



(8) 



(Y 4 5, Y 37 ) = 2(-y y 5 + y x y 4 + y 2 y 7 - y 3 y 6 ) = a u 
(Ym, Y 37 ) = 2(-y 2/6 - 2/12/7 + 2/22/4 + 2/32/5 ) = »13 
(^56, ^37) = 2(2/02/7-2/12/6 + 2/22/5-2/32/4) = a i4 . 



SUB-RIEMANNIAN GEOMETRY OF PARALLELIZABLE SPHERES 



15 



All other vector fields in TL\ U V are orthogonal. For the set Ti 2 we see the following: 

-(Y4 5 ,Y 26 ) = (Y m ,Y 25 ) = -(Y 5e ,Y 24 ) = a 22 
/g\ (^45,^27) = 2(y 2/4 + 2/12/5 - 2/22/6 - 2/32/7) = a 2 i 

(^46)^27) = 2(-y 2/7 + 2/12/6 + 2/22/5 - 2/32/4) = a 24 

(Y56, Y 27 ) = 2(-y y 6 - y x y 7 - y 2 y A - y 3 y 5 ) = a 23 

The other products between vector fields from 7i 2 U V vanish. For 7i 3 the situation is very 
similar. 

(Y4 5 ,Y m ) = -(Y m ,Y 15 ) = (Y 5e ,Y 14 ) = a 33 
/ 10 n (^45,^17) = 2(-y 2/7 - 2/12/6 ~ 2/22/5 ~ 2/32/4) = 034 

' ' (*46, Y X7 ) = 2(-y y 4 + yiyz, - y 2 y 6 + y 3 y 7 ) = a 31 

(Ym, Yn) = 2(-y y 5 - y x y A + y 2 y 7 + y 3 y 6 ) = a 32 

All other vector fields from 7i 3 U V are orthogonal. For the last collection Tt^ we obtain. 

<*45,*0 6 > = -(^46,^05) = (^56,^04) = a 44 

(Ya5,Y 07 ) = 2(y y 6 - 2/12/7 + 2/22/4 - 2/32/5) = «43 
1 ' (Y46, Y 07 ) = 2(-y 2/5 - 2/12/4 - 2/22/7 - 2/32/e) = a 42 

(Y56, Y 07 ) = 2(y 2/4 - 2/12/5 - 2/22/6 + 2/32/7) = «41 
with the rest of the product vanishing. 

We notice some relations between the coefficients CL m k. The coordinates on S possesses 
the equality 

(12) a 2 m + a 2 u + a\ 2 + a\ 3 + a 2 u = 1. 
The direct calculations also show 

«oo + o?i + a\ 2 + a? 3 + a? 4 = 1 

(13) «00 + «21 + a 22 + a 23 + a 24 = 1 

a 00 + a 31 + a 32 + a 33 + a 34 = 1 
2 I 2 I 2 I 2 I 2 "i 
^00 ' ^41 ' *^42 ' ^43 ' ^44 = 

In other words the sum of the squares of the cosines between vector fields from TL m U V, 
m = 1, 2, 3, 4 is equal to 1 — ag . Let us consider 2 cases: < Oq < 1 and a^o = 0. 

Case < oq < 1. 

This case corresponds to any point on S" 4 except of the set 

(14) St = {yl + y\ + y\ + yj = y\ + y\ + y\ + y 2 7 = 1/2}. 
We observe that the sum of the square of the cosines from (I13D : 

4 

^2 a mk = 1 - °oo> m = 1,2,3,4 
fc=l 

belongs to the interval (0, 1) and no one of the cosines can be equal to 1. We conclude that 
each of TC m , m = 1,2, 3,4, is transverse to V. Particularly, if Oqq = 1 then Y^,k=i a mk = 
and H m A-V '. The latter situation occurs in the antipodal points (±1,0,0,0,0) € S 4 or is 
to say on the set 

S2 = {2/o + vl + 2/1 + 2/1 = 0, 2/1 + vl + 2/6 + Vl = i} u 
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(15) {yl + yj + y\ + yl = l, y\ + y\ + yl + y? = 0} e S 7 . 

We also can consider a collection Ho, as a possible horizontal bracket generating distri- 
bution, if we choose an adequate vector field W. If aoo £ (0, 1) we have 

< a 2 A1 + a\ 2 + al 3 + a 2 u = 1 - al < 1 

and no one of the products in (jlip can give 1. We conclude that Yq 7 can n °t be collinear 
to V = {Y±5,Y±q,Ysq}. Therefore, we choose W = Yq 7 . By the same reason we could take 
Yj7, j = 1, 2, 3. In the case when Oq = 1 the vector fields Yj 7 , j = 0, 1, 2, 3 are orthogonal 
to V but Ho is collinear to V and the collection Ho with W = Yj 7 is not transverse to V. 

Case = 0. 

In this case the distribution Ho is nicely serve as a bracket generating if we find a 
suitable vector field W. Notice that (fT2j) became 

(16) ali + a 22 + a| 3 + 044 = 1. 

The a mm can not vanish simultaneously. Without lost of generality, we can assume that 
044 ^ 0. Then a\ x + a| 2 + 043 = 1 — ^44 < 1 from ([13]) and the products (fTTI) imply that 
I07 is transverse to V and can be used as a vector field W. In the case a| 4 = 1 we get 
that Yqj is orthogonal to V. Since WlYj 7 , j = 0, . . . , 3 the collection Ho with any choice 
of Yj7, j = 0, . . . , 3 will be orthonormal. 

We formulate the latter result in the following theorem 

Theorem 4. Let ([5]) be the Hopf map with the vertical space 

V = {Y 45 ,Y 4e ,Y 5e }, 

Si and S2 are given by (fT4"|) and ([7]). Then the Hoph map produces the following Ehres- 
mann Connection H p , p £ S 7 : 

(i) if p $l S\ then H p = (H m ) p , for any choice of m = 1,2,3,4; 

(h) ifp i S 2 then H p = (H U Y j7 ) p , j = 0, 1, 2, 3; 
and we have 

span{7Yo, Yj 7 } p © span{y} p = T p S 7 , j = 0, 1, 2, 3 for all p G S 7 . 

8. Appendix 
8.1. Multiplication of octonions. Let 

ox = (x e + xiei + x 2 e 2 + x 3 e 3 + x 4 e 4 + x 5 e 5 + x 6 e 6 + £ 7 e 7 ) 

and 

02 = (yoeo + 2/1 ei + 2/2e2 + 2/3e3 + 2/4e4 + 2/se5 + 2/6e6 + 2/7e7) 
be two octonions. Then we have according to Table CD 

o\ o o 2 = (a^oeo + £iei + £2e2 + £363 + £464 + £565 + x 6 e 6 + x 7 e 7 )o 
(2/oeo + 2/iei + 2^2 + 2^3 + 2/4e4 + y^ + 2/6e6 + 2/767) = 
= (2:02/0 - £12/1 - £22/2 - £32/3 - £42/4 - X52/5 - x 6 y 6 - x 7 y 7 )e + 
+ (xiyo + xoyi ~ £32/2 + £22/3 - £52/4 + £42/5 + ^72/6 - x 6 y 7 )ei+ 
+ (^22/0 + £32/1 + £02/2 - £12/3 - ^62/4 - £72/5 + £42/6 + x 5 y 7 )e 2 + 
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+(x 3 yo ~ x 2 yi + xiy 2 + x y 3 - x 7 y 4 + x 6 y 5 - x 5 y e + x 4 y 7 )e 3 + 
+(x 4 yo + x 5 y\ + x$y 2 + x 7 y 3 + x y 4 - x\y 5 - x 2 y 6 - x 3 y 7 )e 4 + 
+(xsyo - x 4 yi + x 7 y 2 - x 6 y 3 + x x y 4 + x y 5 + x 3 y 6 - x 2 y 7 )e 5 + 
+ (x6yo ~ x 7 yi - x 4 y 2 + x 5 y 3 + x 2 y 4 - x 3 y 5 + x y 6 + x 1 y 7 )e 6 + 
+(x 7 y + x§yi - x 5 y 2 - x 4 y 3 + x 3 y A + x 2 y 5 - xiy & + x y 7 )e 7 . 

8.2. Right invariant vector fields. According to the previous multiplication rule, we 
have the following unit vector fields of M 8 arising as right invariants vector fields under 
the octonion multiplication. 

Yo(y) = yody + yid yi + y 2 d y2 + y 3 d ya + y 4 d y4 + y 5 d y5 + y e d m + y 7 d yr 
Yi(y) = -yidyo + Vod yi - y 3 d m + y 2 d y . A - y 5 d y4 + y 4 d m - y 7 d ye + y 6 d yr 
Y 2 (y) = -y 2 d yo + y 3 d yi + y d y2 - yid y3 - y 6 d m + y 7 d V5 + j/ 4 ^ 6 - y 5 d y7 
Y 3 (y) = -yzdy - y 2 d yi + yid y2 + y d y3 + 3/7^ + y 6 9 y5 - y 5 d y6 - y 4 d y7 
Y A (y) = -y 4 d m + y 5 d yi + y 6 d m - y 7 d m + y ^ 4 ~ Vidys ~ V2d y6 + y 3 d yr 
Y$(y) = -y<jd m - y 4 d m - y 7 d m - y e d m + y x d y4 + y G d m + y 3 d ye + y 2 d yr 
Y&{y) = -yed yo + y 7 d yi - y A d y2 + y 5 d ys + y 2 d y4 - y 3 d y5 + y d m - yid yr 
Y 7 {y) = -yid yo - y e d m + y 5 d V2 + y 4 d m - y 3 d V4 - y 2 d ys + y x d m + y d yr . 
The vector fields Y{, i = 1,...,7 form an orthonormal frame of T p S 7 , p € S 7 , with 

respect to restriction of the inner product (•, •} from R 8 to the tangent space T P S 7 at each 

pes 7 . 

8.3. Commutators between right invariant vector fields. Denoting by Yij(y) = 
T;[Yi(y),Yj(y)] the commutator between the right invariant vector fields, described in the 
previous subsection, we have the following list: 

Y\ 2 {y) = y 3 d yo + y 2 d yi - yid y2 - y d m + y 7 d y4 + y 6 d y5 - y 5 d y6 - y A d y7 
Yi 3 {y) = -y 2 d yo + y 3 d yi + y d m - y x d yz + y 6 d m - y 7 d y5 - y A d m + y 5 d yr 
Yu{y) = ysd yo + y 4 d yi - y 7 d y2 - y 6 d y3 - y x d y4 - y d m + y 3 d m + y 2 d m 
Yvo(y) = -yid m + y 5 d yi - y 6 d y2 + y 7 d m + y d y4 - yid y5 + y 2 d m - y 3 d y7 
Yie(y) = yid yo + y 6 d m + y 5 d y2 + y 4 d m - y 3 d y4 - y 2 d m - y x d m - y d V7 
Yn(y) = -yedyo + y 7 d yi + y 4 d m - y 5 d m - y 2 d V4 + y 3 d ys + y d y6 - yid y7 
Y 23 {y) = yid yo - y d yi + y 3 d y2 - y 2 d y3 - y 5 d y4 + y 4 d V5 - y 7 d V6 + y e d V7 
Y 24 (y) = y 6 d yo + y 7 d yi + y 4 d m + y 5 d V3 - y 2 d V4 - y 3 d y5 - y d y6 - y x d y7 
Y 25 {y) = -yid yo + y 6 d yi + y 5 d y2 - y 4 d m + y 3 d y4 - y 2 d y5 - y x d y% + y d y7 
Yw{y) = -y 4 d yo - y 5 d yi + y 6 d y2 + y 7 d m + y d y4 + y x d yh - y 2 d ye - y 3 d V7 
Y 27 {y) = y 5 d yo - y 4 d m + y 7 d y2 - y 6 d y3 + yid y4 - y d V5 + y 3 d m - y 2 d y7 
Y 34 (y) = -yjd yo + y 6 d yi - y 5 d y2 + y 4 d m - y 3 d y4 + y 2 d ys - yid y& + y d y7 
Y 3 5(y) = -yedyo - y 7 d yi + y 4 d m + y b d y3 - y 2 d V4 - y 3 d ys + y d m + yid y7 
Y 3 &(y) = y 5 d yo - y 4 d yi - y 7 d y2 + y 6 d y3 + y x d y4 - y d y5 - y 3 d ye> + y 2 d y7 
Y 37 (y) = y 4 d yo + y 5 d yi + y 6 d y2 + y 7 d ys - y d y4 - yid V5 - y 2 d ye - y 3 d V7 
Y 45 (y) = yidyo - y 9 yi - y 3 d m + y 2 d m + y 5 d y4 - y 4 d m - y 7 d m + y 6 d V7 
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^46 (y) = y2d yo + y 3 d yi - y d y2 - ytd y3 + y 6 d m + y 7 d ys - y±d y6 - y 5 d y7 
Y A7 (y) = -y 3 d yo + y 2 d yi - y x d y2 + y d y3 + y 7 d m - y e d y5 + y s d m - y A d y7 
Ym(y) = -Vadyo + V^vx - y\d y2 + y d ys - y 7 d m + y e d y5 - y 5 d m + y 4 d y7 
Y-, 7 {y) = -y2d yo - y 3 d yi + y d y2 + yid ya + y 6 d y4 + y 7 d m - y$ y% - y 5 d V7 

Ye7(y) = yidy - yod yi - y 3 d y2 + y 2 d m - y 5 d m + y A d y5 + y 7 d m - y 6 d y7 . 
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